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ABSTRACT 



We analyze the accretion properties of 21 low mass T Tauri stars using a dataset 
of contemporaneous near ultraviolet (NUV) through optical observations obtained with 
the Hubble Space Telescope Imaging Spectrograph (STIS) and the ground based Small 
and Medium Aperture Research Telescope System (SMARTS), a unique dataset because 
of the nearly simultaneous broad wavelength coverage. Our dataset includes accreting 
T Tauri stars (CTTS) in Taurus, Chamaeleon I, 77 Chamaeleon and the TW Hydra As- 
sociation. For each source we calculate the accretion rate (M) by fitting the NUV and 
optical excesses above the photosphere, produced in the accretion shock, introducing 
multiple accretion components characterized by a range in energy flux (or density) for 
the first time. This treatment is motivated by models of the magnetospheric geometry 
and accretion footprints, which predict that high density, low filling factor accretion 
spots co-exist with low density, high filling factor spots. By fitting the UV and optical 
spectra with multiple accretion components, we can explain excesses which have been 
observed in the near infrared. Comparing our estimates of M to previous estimates, 
we find some discrepancies; however, they may be accounted for when considering as- 
sumptions for the amount of extinction and variability in optical spectra. Therefore, we 
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confirm many previous estimates of the accretion rate. Finally, we measure emission line 
luminosities from the same spectra used for the M estimates, to produce correlations 
between accretion indicators (H/3, Ca II K, C II] and Mg II) and accretion properties 
obtained simultaneously. 

Subject headings: Accretion, accretion disks, Stars: Pre Main Sequence, Ultraviolet: 
Stars, Stars: Chromospheres 



1. Introduction 



Classical T Tauri stars (CTTS) are pre-main sequence objects that are accreting gas from 
their circumstellar disks. The currently accepted paradigm is magnetospheric acc retion, where the 
circumstellar disk is trunc a ted at a few stellar radii by the stellar magnetosphere ( 



I99J iBouvier et al 



1995; 



Johns-Krull et al. 



1999 



Johns-Krul 



^mmm^l ^hhhht — ■ ■ '7 — — — ■ ■ 'J ■ — — ^h^h^h^h^^hh^—^hhh^. 

19981 . 120011 1. with strong magnetic field strengths of a few kG (jJohns-Krull et al 



fc Gafford 



2002 



iartmann et al 



Muzerolle et al 



20001 ). The gas 



at the truncation radius is channeled along the magnetic field lines at nearly free fall velocities, 
until it impacts the stationary photosphere, producing an accretion shock. Accreting sources are 
identified by broad emission lines (like Ha, Ca II and FUV lines of C IV and He II), tracing the fast 



moving material in the accretion streams ([Edwards et al.lll994l ; iMuzerolle et al 



1998 



Ardila et al 



in press! ) . or by an ultraviolet ( UV) and U band excess o ver the stellar photosphere, produced by 
hot gas in the accretion shock (jCalvet fc Gullbringlll998l . hereafter CG98). 



CTTS are known to vary in brightness over short times cales, due to changes in the accretion lu - 
minosity and the modulation of spots on the stellar surface (jHerbst et al.lll994l : lAlencar et al.ll20 12). 



In addition, a spread in t he accretion rate is observ ed for sources of the same age (jHartmann et al 



1998 



Calvet et al. 



2005 



Sicilia-Aguilar et al.ll2010l ). Therefore, the accretion properties of young 



stars are best characterized by studying a large sample, including non-accreting T Tauri stars 
(WTTS) for comparison. WTTS are no longer accreting, but they retain strong magnetic fields, 
due to magnetic dynamo effects, which heat the chromosphere. The ramped up activity i n the 
chromosphere is seen as a UV excess over a dwarf star photosphere (IHoudebine et al.lll996l ) and 
emission in lines of hydrogen and calcium, which are the same d iagnostics of magne tospheric ac- 



cretion; however, the chromospheric contribution is much weaker (jlngleby et al.ll201ll ). Due to the 



similarity in tracers, it is essential to use a WTTS template when estimating the accretion rate onto 
the star (M), especially for low M objects where the accretion emission is comparable to that from 
the active chromosphere. Until now, high signal to noise UV spectra of WTTS were not available 
so dwarf spectra, with lower chromospheric activity levels, took their place as stellar templates, re- 
sulting in overestimated UV excesses. Accurate Ms are vital for the understanding of disk physics. 
The accretion rate provides information about the surface density of the inner circumstellar disk; 
for the lowest accretors, the accretion pro perties reveal the cha r acteristics of the fin al stages of the 
inner disk, shortly before it is dissipated ib'Alessio et"allll999l : Endebv et alJboilh . 
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Initial attempts to fit the UV excess produced in the accretion shock assumed accretion spots 
on the stellar surface wer e characterized by e ither a single tempera t ure and density in treatments 



Rieliaco et al. 


2012 


), or a 


(Calvet & Gullbrine 


1998, 



These models indicated that the accretion spots cover 
a small fraction of the stellar surface, tenths to a few percent. Recent models of the magneto- 
sphere, calculated to reproduce spectro-polarimetric observations, revealed tha t the magnetic field 



geometry of accreting s t ars is complex, incl uding high order and tilted fields (jDonati et ajj|2008j ; 



Gregory &; Donatill201ll ; iGregorv et al.ll2012l ). The accretion footprints on the stellar surface which 
result from field lines of varying strengths and geometries are themselves complex, ranging in size 
and density of the shocked material, with filling factors which may exceed 10% of the stellar surface 
([Long et al.ll201ll ). These large filling factors have yet to be reproduced when fitting the accretion 



excess without significantly over-estimating the UV fluxes. 

In addition to the large filling factors, near infrared (IR) veiling observations are not ex- 
plained by current accretion shock models. Veiling occurs when the excess emission produced 
in the shock fills in, or "veils" , phot ospheric absorption l i nes, causing them to appear shallow 



when compared to a WTTS template (jHartigan et al. 



1999 



literature, out to 1 /im ([Basri fc Batalha 



Edwards et al 



1991 



Valenti et al. 



1993; 



Johns-Krull et al 



Dodin fe Lamzinl I2012T). Veiling at prog r essively redder wave l engths can be found in the 



1990 



Hartigan et al. 



1991 



White &; Hillenbrand 



2006) . More recent results on near IR veiling were discussed in 



Fischer et al. 



2004 



([20 111 ) 



and lMcClure et al.l (lin prepl ). who showed that the amount of veiling becomes nearly constant from 
0.8 to 1.4 fim. For an accretion shock spectrum which peaks in the UV and decreases toward long 
wavelengths, veiling should be negligible in the near-IR. The spectrum of the disk, assuming it is 
produced by circumstellar dust at the dust sublimation temperature around 1400 K, also has little 
contribution to the veiling near 1 fim and does not begin to contribute significantly until ~2 /im 
((Fischer et alihoill ). Long wavelength veili ng may originate in a coo l accretion component, coole r 
than tha t which descr i bes th e UV excess ([Calvet fc Gullbrina Il998l ; IWhite fc Hillenbrand! |2004|) ; 
however, iFischer et al .1 ([201 ll ) found that in some cases the surface area of the cool accretion column 
required to explain the veiling was larger than the stellar surface, so they suggested that the veiling 
came from hot gas in the inner circumstellar disk. On the other hand, by varying the sublimation 
temperature of the dust over a range consistent with different materials, iMcClure et al.l (jin prepl ) 
showed that the veiling between 0.8 and 2.32 microns could be explained by the combined emission 
of dust at the sublimation radius and a cool accretion column with a reasonable filling factor. 

In this paper, we assume that 1 /mi veiling is produced by cool accretion components and use 
the accretion shock models described in CG98 to fit UV and optical observations of a large sample of 
CTTS. The sample is a part of the large Hubble Space Telescope (HST) program, Disks, Accretion 
and Outflows of T Tau stars (DAO; PI G. Herczeg) which compiled a dataset for each source covering 
a long range in wavelength with observations as close to simultaneous as possible. Observations 
included HST far UV (FUV), near UV (NUV) and optical spectra obtained with the Cosmic 
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Origins Spectrograph (COS) and the Space Telescope Imaging Spectrograph (STIS). We focus on 
the NUV and optical observations which are ideal for measuring accretion excesses, whereas the 



FUV spectrum is complicated by molecular e mission, both in lines and the continuum (jBergin et al 



2004 ; llnglebv et al.l l201ll ; France et al.l l201ll ) . Attempts were made to contemporaneously (within 



a few nights of the HST observations) obtain additional optical spectra in order to observe Ha, 
which was not covered in the STIS spectra; these attempts were successful for over half of the 
sample (Section 12 . 3(> . 

Here, we improve upon treatments of the accretion shock as a single spot on the stellar surface 
by including multiple accretion components, a scenario more physically accurate given the geometry 
of the magnetosphere. With the added components, we a ttempt to explain long wavel ength veiling 
by including cooler accretion columns, as suggested by White fc Hillenbrand! ([2004]) and CG98, 
from which the emission peaks toward redder wavelengths. In Section [2] we discuss the data used in 
this paper, both new observations and some previously published. Section [3] describes the process 
we use for calculating accretion rates and our results are presented in Section HI Finally, in Section 
14.11 we explore the range of accretion rates possible, assuming that some accretion luminosity is 
undetectable above the star and in Section [5] we show how our new results affect correlations 
between accretion indicators and M, which are commonly used when UV observations are not 
feasible. 



2. Sample and Observations 
2.1. Description of Sample 



The sample consists of 21 low mass CTTS and 4 WTTS, primarily from the large HST program 
GO 11616 (PI: G. Herczeg), including three CTTS observed earlier with STIS, BP Tau and TW 
Hya (GO program 9081; PI. Calvet) and LkCa 15 (GO pro gram 9374; PI. Bergin). The total DAO 
sample is larger and will be presented in iHerczegl ( in prep ), but here we select those sources which 
are low mass (spectral types late G to early M). The high mass sources in the sample lack accurate 
templates which are necessary for our analysis. The sources inc luded in this paper are in the 
1-2 Myr old Taurus Molecular Cloud (Kenvqn Sz Hartmannlll995l ). the 2-3 Myr ol d Chamaeleon 
I sta r forming region (jGauvin Stroml I1992T), the ^9 Myr 7 7 Chamaeleon region (jLawson et al 



20041 ). the 10 M yr TW Hydra Ass ociation (jWebb et al.lll999l ). and the 16 Myr Lower Centaurus 
Crux subgroup (IPecaut et al.1 120121 ) . Each of these regions is characterized as having isolated low 
mass star formation, with relatively low extinction, so the sources are not affected by remaining 
molecular cloud material or nearby high mass stars. The WTTS and CTTS studied in the paper 
are listed in Tables [Q and [2j 

The majority of the sources are single stars; however there are a few wide binaries. DK 
Tau A, HN Tau A and RW Aur A all ha ve companions at large separations of > 1.4" or ~200 
AU (jKraus et al.1 1201 ll : I White fc Ghezll200ll ). resolved by STIS. CS Cha is a spectroscopic binary 
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(Gucnth er et al.ll20071 ; lNguven et al,ll2012l ). surrounded by a circumbinary disk (jTakami et al.l l2003; 



Nagel et al.ll2012l ). Several of the sources have disks which show evidence for gaps or holes in their 



infrared spectral energy distributio ns, including CS Cha, DM Tau, GM Au r and Lk Ca 15, but all 
of these sources are still accreting (jstrom et al.lll989l : lEspaillat et alJ[2007al lbl l201oh . 



2.2. HST Observations 

Observations were obtained with STIS between 2009 and 2012. STIS NUV observations used 
the MAMA detector and the G230L grating providing spectral coverage from 1570 to 3180 A with 
R~500-1000. Optical observations were completed during the same orbit as the NUV using the 
G430L grating which covers 2900-5700 A with R~530-1040, resulting in almost simultaneous NUV 
to optical coverage with STIS. The bright source, CV Cha, was observed with the echelle grating 
in the NUV (E230M), which we use in place of the G230L spectrum. 

The low-resolution STIS spectra were calibrated with custom written IDL routines following 
the procedures described in the STIS data handbook. The wavelengths were calibrated from the 
location of identified emission lines within the spectrum, and fluxes were calibrated from spectra 
of WD 1337+705 in the NUV and HIP 45880 in the optical. The flux calibration also includes a 
wavelength-dependent aperture correction. Two of th e STIS spectra (those of RECX 1 and RECX 



11) were previously discussed in llnglebv et al.l (|201ll ) 



2.3. Ground Based Observations 

We also obtained low dispersion spectra of much of the sample, covering the Ha line, using 
the Small and Medium Aperture Research Telescope System (SMARTS) 1.5m telescope at the 
Cerro Tololo Inter-American Observatory (CTIO) with spectral coverage between ~5600 - 7000 A. 
Sources were observed with SMARTS for a few days before and after the HST observations when 
possible, providing contemporaneous observations of Ha. The sample of CTTS with SMARTS 
coverage includes AA Tau, CS Cha, CV Cha, DE Tau, DM Tau, DN Tau, FM Tau, PDS 66, RECX 
11, RECX 15, TWA 3a and V836 Tau. For the following analysis, we used the epoch of SMARTS 
spectra whic h was nearest i n tim e to the HST observations. The reduction of SMARTS spectra was 



described in llnglebv et al.l (|201ll ). The flux calibration of the SMARTS spectra was not completely 
accurate due to flux losses in the slit and changes in the observing conditions. To best estimate 
the absolute fluxes, we scaled the SMARTS spectra to the HST optical spectra in the region where 
the two overlap, assuming that the shape of the spectrum is accurate. We assessed the degree 
of variability in the shape of our SMARTS spectra by computing the standard deviation of the 
fluxes in the different epochs and dividing it by the median of the observed spectra. With spectral 
shape variability on the order of 10-25% for our sample (whether due to observing conditions or 
intrinsic variability), there is little error introduced by fitting only one SMARTS optical spectrum 
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per source. 

Due to the low resolution of the SMARTS data, when available, we supplemented the ob- 
servations with non-simultaneous high resolution spectra obtained with MIKE (Magellan Inamori 
Kyocera Echelle) on the Magellan-Clay telescope at Las Campanas Observatory in Chile. MIKE has 
a spectral coverage of 4800-9000 A and resolution of R~35,000. The data were reduced using the 
Image Redu ction and Analysis Facility (IRAF) tasks CCDPROC, APFLATTEN, and DOECSLIT 



(|Todvlll993l ). MIKE spectra were used to calculate V band veiling, which would not be accurate 
at the low resolution of the SMARTS data (Section [372]) . We use MIKE spectra to obtain V band 
veiling for CS Cha, DM Tau, LkCa 15, RECX 11, RECX 15 and TWA 3a. With its location in 
Chile, Taurus observations are difficult and therefore much of our sample was not covered with 
MIKE. Two sources were observed with the echelle spectrograph on the SMARTS 1.5m telescope, 
CV Cha and PDS 66, allowing for veiling estimates. 

Additional observations were obtained with CRIRES, the CRyogenic high-resolution InfraRed 
Echelle Spectrograph on the Very Large Telescope (VLT). CRIRES has a spectral resolution up to 
10 5 in the range of 1 to 5 ^m. These spectra are used in this paper to calculate veiling for a small 
number of sources in our sample with unknown long wavelength veil ing measurem e nts (s ee Section 



3.2f) . A description of the CRIRES data reduction may be found in llngleby et al.1 (120111 ). A log of 
all HST and ground based observations is given in Table El 



3. Calculating Accretion Rates 



3.1. Stellar Template: The Active Chromosphere 

Enha nced chromospheric activity is cha racteristic of stars which have not yet reached the main 
sequence (|Bertoutlll989i ; iGuinan et al.ll2003l ) . A higher level of activity is inferred from lines like Ha 



and t he Ca II infrared triplet, which are in emission and also variable in young stars (jGalvez et al 



20091 ) . WTTS also show exce ss emission at UV through blue optical wavelengths compared to dwarf 
stars (jHoudebine et al.lll996l ). Since Ha, Ca II emission lines and the UV excess are also typical 
tracers of accretion, emission from the active chromosphere must be taken into account when the 
accretion luminosity is estimated from any of the aforementioned indicators. While chromospheric 
e mission is not a large contaminant for strong accretors, it is extremely significant when M is low. 
In llngleby et al.1 (|201ll ) we demonstrated how the emission from an active chromosphere can mask 
all evidence of an accretion shock excess for the lowest accretors, though line profiles observed at 
high resolution revealed complex absorption features produced by the accretion flows. 

Previous estimates of M have main ly relied on dwarf photo spheres as templates against which 



to measure the U band or UV excess (Romaniello et al 



2004). A better estimate of the excess 



could be found using a WTTS template (jValenti et al.lll993l ) but, until now, WTTS with good 
signal in the UV continuum were not available. As part of the DAO sample, WTTS were included 



-7- 



to act as templates, covering the range of spectral types of the accreting sample (Table [Q). Figure 
[1] shows the WTTS templates compared to main sequence dwarf stars with the same spectral type; 
dwarf spectra were scaled to the WTTS at 5500 A. STIS observations of main sequence stars were 
taken from the HST Next Generation Spectral Library (jHeap &: Lindlerl 120071 ). We compared the 
luminosity of the WTTS to that of the dwarf between 2 000 and 3000 A and fo und that the WTTS 
had NUV luminosities ~3x higher than the dwarf stars. iFindeisen et al.1 (|201ll ) showed that WTTS 



of a given age and spectral type exhibit a range in UV fluxes, so the example WTTS shown in 
Figure [T] may not exactly represent the chromosphere for each individual source, but they are the 
best templates to date. In the following analysis, we use the WTTS with the closest spectral type 
match to each CTTS as the stellar template. 
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Fig. 1. — Comparison of WTTS and dwarf stars. In each panel the black line is the WTTS observed 
as part of the DAO sample and the magenta line is a dwarf standard of the same spectral type, 
taken from the STIS Next Generation Spectral Library. The WTTS excess in the NUV is produced 
by an active chromosphere in young stars. 
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3.2. Veiling By Shock Emission 



The excess emission produced in the shock can contribute significantly to the total luminosity of 
the CTTS, making it difficult to determine the luminosity of the star itself. Howeve r, veiled absorp- 



tion l ines provide a diagnostic of the relative contributions from the star and shock (jGullbring et al 



19981 ). Veiling occurs when excess emission is added to the spectrum of the star, filling in the pho 



tospheric absorption lines (IHartigan et al.1 Il989l . Il99ll ) , so by comparing the depth of absorption 



lines in the CTTS to a WTTS, we have an estimate of the veilin g continuum. Th e re are some 
uncertainties in this method; in particular iGahm et al.l (|2008l ) and iDodin fc Lamzinl (|2012l ) show 
that emission lines can also fill in photospheric absorption lines, so treating the veiling emission as 
a continuum may not be correct; however, this effect will likely only be im portant for sources with 



the strongest emission line spectra, or highest veiling (jPetrov et al.l 12011 



The continuum flux responsible for the veiling of absorption lines is added to the flux of 
the photosphere, here taken to be a WTTS, to produce the observed spectrum. So, veiling at 
one wavelength allows us to estimate the intrinsic photospheric flux at that wavelength from the 
observed spectrum. If we have an accurate photospheric template, this scaling provides the stellar 
component flux over the entire spectrum. When available, we used the veiling at V, published in 
Edwards et al.l (|2006l ). which provides a compilation from the literature. For a small sample of the 
sources with unpublished veiling (CS Cha, CV Cha, DM Tau, LkCa 15, PDS 66, RECX 11, RECX 
15 and TWA 3a) we calculated ry from high resolution optical spectra obtained with MIKE or the 
SMARTS echelle. For templates, we used echelle spectra of WTTS with the same spectral type as 
each of our sources. We then added a continuum excess to the WTTS spectrum until the depth of 
the absorption lines matched those in the CTTS, giving us the veiling continuum. We show three 
examples of different degrees of veiling observed in our MIKE spectra in Figure [2j 

Based on the veiling measurements at V band, we scaled our WTTS templates to each CTTS 
using, 

Fv,WTTS = Fv,CTTs/ (1 + r v), (1) 

where Fv,wtts and Fy t cTTS are the continuum fluxes of the WTTS and CTTS at V, respectively, 
and ry = Fy >ve n / Fyy^TTS > is the veiling at V, where Fy jVe u is the excess continuum emission added 
to the photospheric spectrum at V . It is important to note that none of our V band veilings are 
simultaneous to the STIS UV and optical spectra. This introduces uncertainties in our estimates of 
M because CTTS accretion properties are known to be vari able and as M (and the excess continuum 
emission) change, so does the veiling (jAlencar et al.ll2012l ). Ideally, Ms should be measured using 
simultaneous long wavelength veiling information to accurately assess the flux contribution from 
the star and measure the accretion excess. 

According to accretion shock models, veiling should drop to in the near IR as the shock 
spectrum is peaked in the UV (CG98); however, non-zero veiling at long wavelengths has been 
measured in some CTTS, which we attempt to explain in this paper. Veilings a t 1 /im were found 
in the literature for a subset of 10 CTTS in this sample (jEdwards et al.l 120061 ) . For another five 
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Fig. 2. — Veiling in MIKE spectra. The panels show three CTTS (solid, black) observed with MIKE 
which have different degrees of veiling compared to a WTTS of the same spectral type (dashed, 
red). The CTTS are BP Tau (left) with spectral type of K7, CS Cha (middle) with spectral type 
K6 and RECX11 (right) with spectral type of K5. BP Tau has significantly more veiling than 
RECX 11, observed as shallower absorption lines. Typical errors on the veiling are ±0.1. 
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sources, we measured veiling near 1 //m using our CRIRES spectra with either the WTTS RECX 
1 as a template (for CS Cha, LkC a 15 and RECX 11) or using a Near-Infrared Spectrograph 
(NIRSPEC) spectrum of a WTTS (jEdwards et al.l 120061 ). The NIRSPEC template was used to 
calculate 1 fim veiling for DM Tau, FM Tau and IP Tau, after convolving the CRIRES observations 
to the lower resolution of NIRSPEC. For CV Cha, PDS 66, RECX 15 and TWA3a, we did not 
have an appropriat e template to calculate the veiling, so we used the relation, ry = 0.5 x ry 



(Fischer et al 



201: 



L) to obtain ry. We also used this relation to estimate ry for IP Tau, which was 
not in the literature and for which we did we have high resolution optical spectra; however we did 
have CRIRES spectra to calculate ry . All values of ry and ry are listed in Table 



3.3. Extinction and Stellar Parameters 



A large source of error in our estimates of M comes from the assumed amount of extinction, or 
Ay. Our STIS dataset did not provide t he resolution necessary to measure veiling in the optical and 
determine Ay based on the veiling, as in lGullbring et al.l (j 19981 ). so we use extinction estimates from 
the literature. TableHJshows that Ay values for a given source vary widely in the literature. Some of 
the discrepancies in Ay may be due to true variability in the extinction, perhaps by inhomogeneities 
in the intervening molecular cloud or warps in the circumstellar disk. Variable accretion hot spot s 
can also affect t h e colo rs of the star, affecting the Ay determinations (jCarpenter et al.ll200ll . 120021 ). 
Gullbring et al.l (|1998l ) showed that calculating Ay by comparing the colors of reddened stars to 
standard photospheric colors results in different values depending on the colors used. In particular, 
they found that V — I colors are sensitive to the spectral type classification. 

A recently developed method to estimate Ay uses IR veiling estimates to determine the shape 
of the reddened photosphere below the excess and correct accordingly; however, this method has 
only been applied to a few sources in our sample (see Fischer et al. 2011 and McClure et al. 2013). 
Both Fischer et al. (2011) and McClure et al. (2012) note that their estimates are higher than the 
Ays from frequently quoted sources, like lKenvon &: Hartmannl (| 19951 ). There are still discrepancies 
among estimates of Ay for the same source, u sing the IR veiling method (Table HJ). For consistency, 
we use the Ay estimates fromlFurlan et al.l (120091. 120111 ) when correcting Taurus and Chameleon 



I sources for extinction. iFurlan et al 
observed V — I, I — J, or J 



(|2009j, l201ll ). computed extinction values by comparing 



H colors to expected photospheric colors in iKenyon &; Hartmann 



(| 19951 ). using the reddest colors available to minimize the impact from the shock excess which is 
typically strongest in t he blue. For the remaining sources (TW Hya , TWA 3a, RECX 11 and RECX 



15) we assume ^4y=0 (jWebb et al 



1999 



Luhman fe Steeghd 12004 ) . Errors of ±0.5 in Ay can lead 



to an uncertainty up to 1 order of magnitude in M making extinction estimates the largest source 
of error in our calculations. 



The STIS and SMARTS spectra were de-reddened using the reddening law of IWhittet et al 



(|2004l ). UV emission is extremely sensitive to reddening assumptions, with A\ > 2 x Ay nea r 
2500 A and corrections are complicated by uncertain UV reddening laws. ICalvet et al.l ( 2004 ) 
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compared UV reddening laws and found that of lWhittet et al.l ([20041 ) was appropriate for sources 
in environments like the Taurus Molecular Cloud. We calculated the s tellar luminosity from t he flu x 
in the J band of the photosphere using the bolometric correction of iKenvon fc Hartmannl ([1995 ). 
The photospheric J band flux was obtained by scaling de-reddened 2MASS J magnitudes by the 
veiling at 1 fim. Although there is an excess above the photosphere at J, with known veiling we 
separate the J flux of the star from that of any continu um excess and use the flux from the star 
alone when calculating the luminosity. We then used the ISiess et al.l (J200CJ) evolutionary tracks to 
determine the masses of the sources in the sample and estimate radii from the luminosities. The 
stellar properties for our WTTS and CTTS are listed in Tables Q] and [2J respectively. 



3.4. Accretion Shock Model 

A full description of the model used to characterize the accretion shock can be found in CG98 
but here we review the main points. The current picture of material accretion in the inner disk 
of CTTS is magnetospheric accretion. Columns of accreting material fall onto the star along the 
magnetic field lines traveling at the free fall velocity, v s , hit the stationary photosphere and create 
a shock. The velocity of the material is given by 

where M* and R* are the stellar mass and radius, respectively, and Ri (assumed to be 5 R*) is the 
radius at which the magnetosphere truncates the disk (CG98). 

The model is simplified by assuming that the accretion column has a plane parallel geometry 
and is perpendicular to the stellar photosphere. The shock formed at the base of the accretion 
column reduces the velocity of the infalling material in order for it to join the star at the photosphere, 
converting the kinetic energy into thermal energy, and causing the temperature to increase sharply. 
In the strong shock approximation, which is used because the material is traveling at high velocities, 
the temperature immediately after the shock is described by 

(CG98). At these high temperatures, the shocked material emits soft X-rays into the pre-shock, 
post-shock and the photosphere below the shock. This radiation heats the material in these regions 
causing it to emit the observed excess continuum emission. 

In this treatment, the emission from the accretion column is characterized by two parameters, 
T and /, which are the total energy flux in the accretion column [T = l/2pVg, where p is the 
density of the material in the accretion column) and the filling factor (or fraction of the surface 
of the star which is covered by the accretion hot spot), respectively. Changes in the filling factor, 
/, cause the shock spectrum to increase or decrease independent of wavelength, in essence scaling 
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the luminosity of the emission. T acts to change the wavelength of the peak of the accretion 
shock emission. Figure [3] shows how the shock spectrum, normalized by the maximum flux of each 
spectrum, changes for different values of J 7 , assuming typical values for M* and i?* of 0.5 M & and 
2.0 Rq, respectively. As T increases, more energy is deposited on the stellar surface, increasing 
the temperature of the photosphere below the shock. This region behaves as the photosphere of a 
star with an earlier spectral type than that of the undisturbed photosphere. High T columns have 
temperatures up to 9000 K, whereas low T columns are cooler, around 5000-6000 K, producing 
the wavelength shift between columns with different J-s. For the lowest value of J-, the emission 
peaks around 6000-7000 A whereas for the model with the highest the emission peaks between 
2000 and 3000 A. 
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Fig. 3. — Fluxes of accretion shock models with varying energy flux {T). Each line shows the 
emission from an accretion column with a different value of J-, normalized to unity. The energy flux 
of the column is listed next to each spectrum, spanning the range from T = 10 10 — 10 12 erg s _1 cm -2 . 
The peak of the emission shifts to longer wavelengths as T decreases. 

We use this property of the accretion column emission to explain excesses at both short and 
long wavelengths when attempting to fit the DAO observations. In the following text we refer to 
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accretion columns with energy fluxes of T < 10 11 erg s _1 cm -2 as "low J 7 " columns and those 
characterized by T > 10 11 erg s _1 cm~ 2 as "high J 7 " columns. J- depends on the density and 
velocity of the accreting material; we assume that the velocity is constant, as set by the geometry 
(see Equation [2]) , so the range of T represents high and low density columns. Note that as T 
decreases, larger values of / are needed to produce the same amount of flux as a higher T accretion 
column. M can be calculated for each accretion column, with T and / known, using 



M 



8irRi 



(4) 



(jGullbring et al.l l2000l ) . The total M is a sum of the contributions from each column and from M 
we determine L acc using, 



GA'LM 



R* 



:i 



(5) 



This treatment does not distinguish between multiple accretion spots, each with a distinct density, 
or a single accretion spot with a range of densities. The latt er scenario is suggested by models pro- 



duce d using the Zeeman Doppler Imaging technique (ZDI) ([Gregory &: Donatill201ll : iDonati et al 
20081 ). 



In this paper, the excess which veils the photospheric absorption lines is primarily a contin- 
uum. The main opacity sources included in the calculation are H bound- f ree and free- f ree, H - 
bound-free and free-free, C, Si and Mg, plus additional sources describ e d in Calvet et al.l (| 19911 ). 
Line blanketing is also included, using the line list from Eciirucz &; Belli ( 19951 ). A more thorough 
treatment of the s pectr um of the excess continuum including spectral lines has been presented in 
Dodin Sz Lamzinl (|2012j) where it was found that emission lines are important to consider when 
estimating veiling. Dodin & Lamzin showed that veiling may be overestimated when lines are 
neglected in the veiling calculation. An important difference between the models of CG98 and 
Dodin &: Lamzinl f|2012l ) is the structure of the region where the heated photosphere joins the shock 
above it, which in both treatments is taken at the ram pressure of the shock. In CG98, when the 
heated photosphere is joined to the postshock region the temperature rises quickly to ~ 10 6 K 
because this region has small physical dimensions. Metals would become quickly ionized and few 
photospheric emission cores would be expected, especially s ince the location of this transition re- 
gion is close to the temperature minimum. In contrast, since iDodin fe Lamzinl (120121 ) do not have a 
rapid temperature rise, strong emission lines form in this region. Additional work is needed to sort 
out the differences between the two models. For now, we acknowledge that emissi on lines may affect 
our results, especi ally for the highly veiled sources like DR Tau and RW Aur A (jGahm et ajJ|2Q08j ; 
Petrov et allboilh . and accretion spot sizes for those sources, in particular, may be overestimated 
due to the omission of emission lines in the veiling spectrum. 
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4. Results 

In this analysis we include, for the first time, multiple accretion columns covering a range in 
energy flux; specifically, we add low T accretion columns to high J- columns to fit both the UV 
and optical excesse s. Evidence for multiple accretion spots is seen in maps of the magnetosphere 



(jDonati et al.ll2008r) and also in fits of the broad He and Ca emission, which require inhomogeneity 



in the hot spots (IDodin et al.ll2013l ). We also compare the amount of 1 \xm. veiling predicted by 
our models with observed values. Our spectra do not extend into the near IR, so to estimate 
the 1 fim veiling produced by our model we first find J of our WTTS template. We scaled V 
by the veiling, ry = Fy tVe u/ 'Fv^wtts, and then using V — J colors for the given spectral type, 
find the J magnitude of the WTTS which reflects the veiling. We assume tha t r.j = n jim , which 



is valid because the v eiling continuum is constant between 0.8 and 1.4 /im ([Fischer et al.l 12011 



McClure et alJIin prepl ) . We estimate J of the model by calculating t he flux that would be measured 



assuming the transmission curve for the Two Micron All Sky Survey (jSkrutskie et al.ll2006l . 2MASS) 



To find the best fit of the accretion shock models to the STIS and SMARTS spectra for each 
source, we calculated the emission from accretion columns spanning T = 10 10 - 10 12 erg s 1 cm 3 , 
using the masses and radii from Table [2J We allowed the filling factor, /, to vary independently 
for each accretion column and then summed the flux from each column with that of the WTTS, 
scaled to the CTTS by the veiling at V (which tells us the location of the photosphere below the 
excess). The spectrum of the final model is given by; 



\,model 



Y,i(F x (column,Ti) x /(J^)) + f x,wtts, (6) 



where F\ mo dei is the flux in the model, F\(column, Fi) is the flux from an accretion column with 
energy flux T{ and f{J-%) is the filling factor of a column with energy flux Ti, where J^=10 , 3 x 
10 10 , 10 , 3 x 10 , 10 12 erg s _1 cm -3 . The values of T we have chosen represent the typical 
range of values for CTTS (CG98). 

When determining the best fit model, we isolated continuum regions because the accretion 
shock models do not attempt to reproduce line emission produced in the shock. We calculated the 

2 

Xred' as 

2 _ 1 V JV {F\,model ~ F\^ C TTs) fff . 

Xred — J7^i=0 vJ 

where is the number of continuum wavelengths which contribute to the fit and E\^ctts is the 
error in the observed fluxes. The filling factor of each contributing column in the best fit model is 
given in Table [5j columns with / = do not contribute to the model fluxes. Finally, we calculated 
M using Equations [21 H] and, 

M^xS^x/^)) (8) 
v 

Estimated Ms are listed in Table along with the sum of the filling factors of all columns. 
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The best fit of the accretion shock models to the de-reddened NUV and optical spectra, cor- 
responding to the stellar and accretion properties listed in Tables [2] and [5l are shown in Figures 
H]-[8j Overall, we found good fits between the models and the observed spectra, with a few ex- 
ceptions. For the later spectral types, especially the M stars, the models did not reproduce a rise 
in the observed spectra between 2000 and 3000 A (see, for example, RECX 15 and DE Tau in 
Figure E]). This spectral region is populated with Fe emission lines, which a re unresolved at the 



STIS resolution. T he Fe lines may be produced in accretion related processes (jHerczeg et al 



2005 



Petrov et allboilh : however, they are also observed in the WTTS templates, though weaker, they 



appear to have a chromospheric component. Fe lines become more apparent at later spectral types, 
as the photospheric emission in the UV decreases. Finally, we do not attempt to fit the FUV 



20091 ) that are not 



spectrum (<1700 A) which has contributions from H2 in the disk ( Inglebv et al. 
included in the accretion shock model. 

Figures 0]- [8] also include non-simultaneous photometry from the literature, de-re ddened usinf 
the Ay s from Table [2j When available, we used the range of optical photometry from lHerbst et al.l 
(|1994l ) and this range is shown as the gree n error bars. The p hotometry for the remaining CTTS 
came from the follo wing sources; CV Cha (ILawson et al.lll996T). GM Aur, IP Tau, LkCa 15, TW 



Hya and V836 Tau (IKenvon Ha rtmann 1995 ), PC' S 66 (IBatalha et al 



RECX 11 and RECX 15 (|Sicilia-Aguilar et al 



1998 



Cortes et al.l 



20091 ) and TWA 3a (|Gregorio-Hetem et al 



2009) 



1222) • 



The models shown do not attempt to fit the photometry because it is not simultaneous, but we use 
the photometry to look for evidence of high amplitude variability. Most of the photometry agrees 
with the STIS and SMARTS spectra and is therefore well fit by the models. For a number of the 
sources, RW Aur A, FM Tau, IP Tau, V836 Tau, RECX 15 and TWA 3A, the STIS spectra are 
slightly lower than expected from the photometry. Of these sources with mis-matched photometry 
and spectra, with the exceptions of RW Aur A and FM Tau, photometry for only one epoch is 
available, so the difference may be due to variability. The DAO N UV and optical spe ctra of HN 
Tau A are significan tly higher than the range of photometry from iHerbst et al.l (|1994l ); however, 
Grankin et al.1 (|2007l ) observed that the brightne ss of HN Tau A was steadily decreasing during the 
10 years between 1985 and 1995. It may be that lHerbst et al.1 ()1994l ) observed HN Tau A during a 
dim period and the source has since brightened; su ch changes in brightness are observed in other 
CTTS, so it is not unexpected (jGrankin et al.1 120071 ). 



In Table [5] we also give the values of ry predicted by our models. Comparing the model ry to 
those from the literature or our observations (Table [2|), the model values for all but four sources 
are within ±0.1 of the observed ry, indicating that near IR veiling naturally occurs from a multi- 
column accretion shock model. Sources with large ry are fit by models with large filling factors 
for the low T columns. The four sources which do not agree are AA Tau, DR Tau, DK Tau A 
and HN Tau A. For DK Tau A, we were unable to reproduce the 1 ^m veiling, likely because of 
variability. In lEdwards et al.1 (|20061 ). the veiling at V is equal to the veiling at 1 /im, requiring an 
abnormally red accretion spectrum, which is not supported by our NUV and optical observations. 
The V and 1 /im veiling measurements were not simultaneous, so it is possible that M was higher 
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when the IR observations were obtained. HN Tau A also appears to be a remarkably variable 
object, as mentioned above, so it is possible that variability is responsible for not being able to fit 
the non-simultaneous optical and near IR veiling data. Veiling is very difficult to measure for DR 
Tau because the excess continuum is so strong that many absorption lines are completely filled in; 
therefore, we expect large error bars on the near IR veiling. 

The models in Figures 0]- [8] represent the best fit to the UV and optical data; however there is 
some degeneracy in J- and /. There is a range in T and / which will produce a shock spectrum that 
results in a comparable x 2 re d ^ *° the data. Figure [9] shows this range for CS Cha and FM Tau, 
where the colored points represent values of Xred ^ or m °dels with different total filling factors and 
characteristic J- values. The characteristic J- value is the average T of the contributing columns, 
weighted by the filling factor of each column. Figure [9] reveals that there is a range of roughly one 
order of magnitude in both J- and / for which models will result in fits with Xred < 2x the best 
fit model. While interesting for the shock geometry, this degeneracy has little affect on the M's 
because each model is accurately measuring the total excess from which the accretion luminosity 
is measured. 

An important result is that, by including the low T columns, we find significantly higher filling 
factors than would be needed to fit the UV excess alone (Table [5]). CG98 fit optical spectra of 
Taurus CTTS using a single accretion column model. For seven out of nine of our sources which 
overlap with CG98, we find higher / values. The characteristic T of the contributing accretion 
columns is typically lower than that found by CG98, because of the addition of low T accretion 
columns which have large filling factors. However, both DN Tau and AA Tau have larger filling 
factors and higher T values in CG98 than in our analysis. DN Tau was fit by an accretion column 
characterized by log J- = 10.5 and / = 0.5% in CG98; however, their fit to the observed spectrum 
cannot explain the blue excess very well. Our model includes a high T column along with the 
low T column, better fitting both the blue and red excesses and skewing the characteristic T 
to a higher value. We mentioned above that our model for AA Tau did not reproduce near IR 
veiling measurements. There may be variability in the cool accretion component, where it was less 
prevalent during the DAO observations. Were we to increase the contribution from low T columns 
to fit the 1 /um veiling, it would bring the characteristic J 7 to a lower value and the f to a larger 



value than found in CG98. Given that models of the magnetosphere predict large / (jLong et al 



2011 



), multi-component accretion columns are a better representation of the physical geometry of 
the system and now reproduce the flux distribution as well. 

We compare the accretion properties calculat ed here with those fr om two previous studies of 



young stars in Taurus in Figure [10] and Table [U IValenti et al.l (119931 ). hereafter V93, fit optical 
spectra with coverage from 3400 to 5000 A of Taurus CTTS using a hydrogen slab model and a 
WTTS template. The slab wa s characterized by the temperature, number density, thickness and 



coverage on the stellar surface. iGullbring et al.1 (119981 ) . hereafter G98, used similar assumptions to 



those of V93 and therefore, their accretion rates are in good agreement. G98 notes that the major 
differences between their Ms and the V93 Ms come from choices of Ay and the evolutionary tracks 
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Fig. 4. — Spectra of late G and early K spectral type CTTS in the DAO sample: CV Cha, PDS 66 
and RW Aur A. We fit the STIS NUV and optical spectra with emission from accretion columns 
plus a WTTS spectrum. In each panel the black spectrum is the CTTS and the blue spectrum is the 
WTTS, LkCa 19. The broken black lines represent accretion shock models with different J- values, 
defined as in Figure [3l The red line is the best model fit to the data, adding the emission from the 
different shock models to the WTTS spectrum. The green error bars indicate non-simultaneous 
photometry, representing the range observed in multi-epoch observations when available. Although 
PDS 66 does not have an NUV excess, examination of the Ha emission line reveals that it is broad, 
and therefore accretion is likely ongoing, though at low levels. 
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Fig. 5.— Spectra of mid K CTTS in the DAO sample: CS Cha, DR Tau, HN Tau A, LkCa 15 and 
RECX 11. Lines are defined as in Figure H The adopted WTT S in each pane l is RE CX 1. RECX 
11 does not have a significant NUV excess, but it was shown in llnglebv et al.l (|201ll ) that it is still 
accreting. 
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Fig. 6. — Spectra of late K and early M CTTS in the DAO sample: AA Tau, BP Tau, DK Tau A, 
DN Tau. Lines are denned as in Figure HI The adopted WTTS in each panel is HBC 427. 
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Fig. 7. — Spectra of late K and early M CTTS in the DAO sample continued: FM Tau, GM Aur, 
IP Tau, TW Hya and V836 Tau. Lines are denned as in Figure HI The adopted WTTS in each 
panel is HBC 427. 
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Fig. 8.— Spectra of mid M CTTS in the DAO sample: DE Tau, DM Tau, RECX 15 and TWA 3a. 
Lines are denned as in Figure [U The adopted WTTS in each panel is TWA 7. 
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Fig. 9. — Degeneracy in characteristic J- and total filling factor f. The color of each point represents 
the value of x 2 re d calculated from Equation assuming a model with a given total filling factor and 
characteristic J- (the average J- of the contributing columns weighted by the filling factor) for CS 
Cha (left) and FM Tau (right). The plot range represents a typical parameter space over which the 
accretion shock models were calculated. The red box shows the point where Xred 1S a ^ a minimum. 
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used to determine M*. Different Ays change the stellar luminosity, and therefore radius, so the 
value of M^jR* will be affected by the extinction estimate. To avoid the additional variables of 
mass and radius in the analysis, we compare the accretion luminosities from each study instead of 
M. 

In Figure HUJ left, we plot our values of L acc compared to those of G98 and find that our 
estimates tend to be higher. The values of L acc we calculate here (shown by the asterisks) have less 
than a 30% difference from those of the G98 for half of the sample. After correcting our accretion 
luminosities by the difference in Ay, shown by the arrows, we find that all but one are in good 
agreement. The outlier is HN Tau A, where our M is over one magnitude larger than G98. This 
is likely due to real variability; we observed that t he STIS fluxes are higher than the range in 
optical photometry observed by iHerbst et al.1 ()1994l ). It appears that two factors in our analysis, 
the excess flux decrement in the UV due to the WTTS template and the additional excess flux in 
the red due to the low J- column, cancel each other out. Therefore, we find similar values of L acc , 
af ter accounting for d i fferen ces in Ay, to those of G98. We chose to use the extinction estimates 
l|200flL 



of Fur lan et al. 



2011 



who used IR colors to esti mate the amount of exti nction, avoiding 



spectral regions where the shock contribution is highest. iFurlan et al. (boodhonh also covered all 
of our sources in the Taurus and Chamaeleon I regions, allowing us to use Ay's which were derived 
consistently for each source. 

When comparing the DAO accretion luminosities to those of V93, our estimates are high for 
three sources; DK Tau A, HN Tau A and DN Tau and low for two; FM Tau and RW Aur A. V93 
obtained optical spectra from the UV Schmidt spectrometer at Lick Observatory and provided an 
atlas of the observations making direct comparison of the optical fluxes possible. For DK Tau A, 
HN Tau A and DN Tau, the optical fluxes are lower in V93 than in our observations while for RW 
Aur A and FM Tau they are higher, consistent with the discrepancies in L acc . These variations in 
the observed fluxes may be intrinsic variability but the difficulty of flux calibrating ground based 
slit spectra may also contribute. In particular for FM Tau, V93 notes that their slit loss correction 
was greater than 25%, indicating that it was observed at a high zenith angle or in poor seeing, 
making flux calibration more uncertain. 



4.1. Hidden Accretion Emission 



Not all sources require a low energy column to fit the UV and optical excesses in our analysis; 
however, there could be emission from cold columns hidden by the photospheric flux. If these 
columns were present, they would increase our estimates of the mass accretion rates. Here, we 
determine how much hidden flux may be present for sources with no detectable red excess. In 
Figure [TTI we show two mode l fits to the spectrum of V836 Tau, which has no observed veiling 
at 1 fim (jEdwards et al.ll2006l ). In the left panel, the best fit model is produced with only a high 
T column, no additional accretion columns are needed (Table [5]). In the right panel, we assume 
that a low T column exists, but the emission from the column is not detectable above the stellar 
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Fig. 10. — DAO accret ion luminosities v ersus values from the lit erature. We c o mpar e our estimates 
of L acc to those from Gullbring et al. ( 1998 ) on the left and Valenti et al. (jl993l ) on the right. 
The solid line shows where sources would fall if the accretion luminosities agreed. The dashed 
lines represent a ±30% difference in the values of L acc . Asterisks represent accretion luminosities 
calculated assuming th e value of Ay listed in T able [2] while the arrows show how L acc would change 
if we use the Ays from I Gullbring et al.l (| 19981 ) or lValenti et all (j 19931 ) . 
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component. For V836 Tau we find that the contribution to the accretion emission that may be 
hidden below the stellar emission could be equal to that in the high T column, doubling the 
estimated accretion rate. This new model meets the constraint of having limited veiling at 1 /im, 
with n < 0.1, within the errors of veiling estimates. 

We perform this analysis for all sources which had / = for the low J- columns in the initial 
fit, with the constraints that the NUV and optical fluxes are not overestimated by the new models, 
and that the veiling at 1 /iin remains within ±0.1 of the observed 1 /im veiling. Using this method, 
we calculate upper limits on M and the filling factor, listed in Table [TJ where the differences in the 
two filling factors comes from increasing /(l x 10 10 ) from its value in Table [5j These additional 
low J- accretion columns are important because they increase the expected M along with the filling 
factors. For the example shown in Figure \TT\ the accretion rate would double if we include the 
maximum shock emission that could be hidden by the photosphere, but the filling factor may be 
up to 100 times higher. The increase of the filling factor is much higher than M because there is 
less mass per unit area in the low T columns than in the high T columns. This would increase the 
filling factor from <0.1% to a few percent, better in line with the m odels of the accretion footprints 
produced by complex magnetic field geometries (jLong et al.ll201ll ) and the distribution of excess 



accre tion related emission regions now being found from magnetic mapping studies (jDonati et al 



2011 



) . The addition of a hidden accretion column is increasingly significant for sources with a lower 
accretion rate, where the flux that may be hidden by the star becomes comparable to the excess 
observed in the UV. 
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Fig. 11. — Range of possible accretion shock models for V836 Tau. In each panel the lines are 
defined as in Figure HI The left panel shows the best fit model of the accretion shock to the optical 
spectra. The right panel shows the fit if we assume that some accretion flux is hidden by the 
stellar photosphere. The upper limit on the luminosity of the low T accretion column (dotted 
line), constrained by the lack of veiling at 1 p, is equal to that of the column which fits the NUV 
excess (long dashed line), doubling M. 
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Correlations with Accretion Indicators 



Measuring the UV excess is ideal for calculating all but the lowest Ms; however, UV observa- 
tions are difficult to obtain. For this reason it is common to use an indicator of accretion easily 
accessible in the optical or near-IR which has been calibrated with Ms calculated for the few sources 
with UV observations. We assume that broad optical line emission, primarily in lines of hydrogen 
and calcium, originates in the material free-falling onto the star in the accretion flows supported 
by the fact that the emissio n lines are well reproduced by models which assume this geometry 
(|Muzerolle et al.lll99sl . hoOlh . Recently, iDupree et al.1 (|2012l ) sugg ested that the lines are instead 
formed in the turbulent region directly below the shock. However, ICalvetl ()198ll ) found that even 
a chromosphere covering the entire stellar surface did not have enough emitting volume to explain 
the lower Balmer lines in CTTS; the shock region also has a small emitting volume and so it is 
unlikely to explain the observed emission in these lines. 

Both the equivalent width of Ha, EW(Ha), and the width of Ha at 10% of the maximum 
flux are often used as accretion proxies (I White &: Basrill2003l : iNatta et al.l 12004 ). The width of Ha 
at 10% is expected to be a better indi cator of accretion, sin ce the wings of the line trace the fast 
moving material in the accretion flows (jMuzerolle et al.ll200lh. Also, the E W(Ha) saturates for the 
highest accretors, or sources with the largest veiling (jMuzerolle et al.lll998l ). Correlations between 
M and Ha 10% width have been widely used; however there is considerable scatter because Ha 
and the UV excess have not been measured simultaneously. The dataset presented here, with UV 
and Ha observations separated by less than one to a few days, would be ideal for this analysis; 
however, the resolution of the SMARTS optical spectra is too low to accurately measure the 10% 
width of Ha. Even ob serving the wings of Ha may not be the best method for determining M, as 
Costigan et al.1 (120121 ) showed that variability in the wings exceeded variability observed in other 
tracers. 

In Figure \12\ M and EW(Ha) are plotted for the sample of our sources which had nearly simul- 
taneous SMARTS optical and UV observations. There is no clear relation between EW(Ha) and 
M. Two sources, DM Tau and RECX 15, have relatively low M but sign ifica nt EW(Ha), showing 



that it is not a good tracer of M. However, iHerczeg &: Hillenbran d (2008) and 



Manara et al 



Herczeg fc Hillenbrand! (|2008h 



(2012) 



found correlations between the luminosity of Ha and M or L acc 
determine L ar a using optical spectra with coverage down to the atmospheric limit in the blue while 
Manara et al.1 (|2012l ) observed the U excess (to estimate L acc ) and Ha with the HST Wide-Field 



Planetary Camera 2 (WFPC2). Our analysis is the next step in defining these correlations because 
we have long wavelength spectral coverage, in particular the STIS spectra provide the shape of the 
excess in the UV. We can measure the flux of the Ha line in the SMARTS data and Lh« is then 
calculated using distances from Table [2j In Figure [13] we see a correlation between Ha luminosity 
and M (or L acc ) with a Pearson correlation coefficient of 0.9. The least-square fit to the data in 
Figure [13] yields 

log(M) = l.l(±0.3)log(L Ha ) - 5.5(±0.8) (9) 
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log(L acc ) = 1.0(±0.2)log(L Ha ) + 1.3(±0.7) (10) 

For Equations [9Hl6| M is in units of M yr -1 and all luminosities a re in units of L & . Equation [10] 
is in good ag reement with t he fit s found in lHerczeg &: Hillenbrand! (j2008l ) and differs only slightly 
from that of iManara et al.l (|2012l ). who found a y-intercept of 2.6±0.1 as opposed to our value of 
1.3±0.7. 
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Fig. 12. — M versus EW(Ha). We compare our calculated Ms given in Table [5] to the EW(Ha) 
measured from the nearly simultaneous SMARTS spectra. The asterisks show the initial fits to the 
UV and optical excesses. Arrows show the range of M for sources which were assumed to have a 
hidden cool accretion component, with the point of the arrow representing the upper limit on M. 
We see no correlation between the EW(Ha) and M. 



H/j is less often used to trace accretion than Ha but has bee n shown to correlate with M 
dMuzerolle et alJboOll : iHerczeg & Hillenbrandlbood : iFang et alibood ) . We compare the luminosities 
in H/3 (Lh/3) measured from the simultaneous STIS optical spectra to both M and L acc in Figure HH 
We find strong correlations between the line luminosity and the accretion rates and luminosities, 
both with correlation coefficients ~0.9. The lines shown in Figure Q3] which describe the trends are 
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log L(Ho, L @ ) log L(Ho, L ) 

Fig. 13. — M and L acc versus Ha line luminosity. The luminosity of Ha, measured in the low 
resolution SMARTS spectra, is correlated with the accretion properties. Asterisks and arrows are 
defined as in Figure [12j Least square fits to the data are given in Equations [9] and [10l 
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given by the following equations; 

log(M) = 0.9(±0.1)log(L H/3 ) - 5.1(±0.5) 
log(L acc ) = 1.0(±0.1)log(L H/3 ) + 2.4(±0.5). 



(11) 
(12) 



Again, our relation between L acc and Lh@ agrees with that found by lHerczeg fc Hillenbrand! (j2008l ). 




Fig. 14. — M and L acc versus H/3 line luminosity. Asterisks and arrows are defined as in Figure 
[T21 We measured the luminosity of H/3 from the simultaneous STIS spectra. H/3 clearly traces the 
accretion properties of the source. Fits to the data are given in Equations [TT] and [T21 



Another commonly u sed tracer of accretion is the Ca II infrared triplet line emission (jMohantv et al 



20051 : iRigliaco et al.ll201ll ). Our SMARTS spectra did not cover the infrared triplet lines of Ca II, 
but in our STIS optical data, we observe the Ca II H and K lines. The H line at A3969 A is 
blended with the A3970 A He line at the resolution of STIS, so we only consider the K line which 
is fr ee of contamination. Both the H and K lines are also used as tracers of chromospheric activ- 
ity (jMamaiek &: Hillenbrand! 120081 . and references therein). We measured the luminosity of the K 
line in our sample of WTTS and found that the maximum value of log L(Ca II K) = —4.3 L@. 
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For CTTS with similar Ca II K luminosities, there will be some contamination in the line due to 
chromospheric activity. Similar to the Ha and H/3 line luminosities, the Ca II K line at A3934 A 
appears to have an origin in accretion related processes, as there is a clear trend between the Ca II 
K luminosity (Lcaii k) and M or Lace (with, correlation coefficients of 0.8 in each. case). We show 
the trends in Figure [T5l and give the relations in the following; 

log(M) = 0.9(±0.2)log(L C aii k) - 5.1(±0.7) (13) 
log(L acc ) = 1.0(±0.2)log(L C aii k) + 2.2(±0.7) (14) 
E quation [T4"l is once again i n goo d agreement with the relation between L acc and Lc a ii k produced 



m 



Herczeg fc Hillenbrand! (|2008l ). 
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Fig. 15. — M and L acc versus the Ca II K line luminosity. Asterisks and arrows are defined as in 
Figure [12j We measured the luminosity of the Ca II K line from the simultaneous STIS spectra. 
Like the infrared triplet lines of Ca II, the K line luminosity correlates with the accretion properties. 
Fits to the data are given in Equations [13] and Q3J 



Included in our NUV spectra are tw o emission lines which have been shown to correlate with 
accretion properties. ICalvet et al.l (|2004l ) analyzed NUV spectra for a sample of intermediate mass 
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T Tauri stars and found that luminosities of C II] A2325 A and Mg II A2800 A increased with L acc . 
We expected a similar trend due to simultaneous observations of the lines and accretion excess and 
Figure [161 shows a clear correlation. Mg II is also a tracer of chromospheric activity, and is detected 
in each of the WTTS templates, but the chromospheric contrib ution to the luminosity is likely 
saturated at the age of our sample (|Cardini fc Cassatellal 120071 ) and therefore not contributing 
to the spread in luminosity. In the CTTS, Mg II can have strong absorption features due to 
the presence of outflows, which was seen in a sample of NUV spectra of T Tauri stars obtained 
with the International Ultraviolet Explorer (jArdila et al.l |20Q2| ) . C II] is not observed in WTTS; 
however, it is easily identified even in low M objects and correlates with L arr . Therefore, C II 



appears to be a very sensitive tracer of accretion or accretion related outflows (jCalvet et al 



2004 



Gomez de Castro &: Ferro-Fontanl 120051 ) . observable even when the NUV continuum excess is not. 
Relations between line fluxes and L acc are given by; 



log(L acc ) = l.l(±0.2)log(L c n] ) + 2.7(±0.7) 
log(L acc ) = l.l(±0.2)log(L Mg n) + 2.0(±0.5) 



(15) 
(16) 



In conclusion, there are several secondary tracers of accretion which are easily accessible in 
optical spectra and the good agreement between the correlations found here and in the literature, 
using different techniques, show that these correlations are robust. With well calibrated accretion 
indicators, it is not necessary to go to the UV to obtain accretion rate estimates. These results will 
facilitate the study of accretion properties for large samples of sources using optical spectra, which 
are significantly easier to obtain than UV data, and may offer more opportunities to measure 
M simultaneously to data tracing other T Tauri phenomena, like outflows or circumstellar disk 
dynamics. 



6. Summary and Conclusions 

We used a large sample of T Tauri stars with nearly simultaneous spectral coverage in the 
FUV, NUV and optical to analyze the accretion properties of young stars. The main results are 
summarized here. 

• The NUV fluxes of WTTS are ~3x stronger than dwarf stars of the same spectral type, 
due to enhanced chromospheric activity. This emission can introduce uncertainties in the 
estimate of the NUV excess produced by accretion if the correct templates are not used. By 
using WTTS templates with good signal to noise in the UV, which first became available 
in the DAO sample, we are able to distinguish chromospheric excesses from accretion excess 
when calculating accretion rates. 

• Accretion shock models which explain the NUV excess cannot describe both large filling 
factors predicted by models of the magnetic field distribution and non-zero 1 /im veiling. 
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Fig. 16. — Accretion luminosity versus luminosity of NUV emission lines. We compare the lumi- 
nosity of C II] A2325 A (left panel) and Mg II A2800 A (right panel) to the range of accretion 
luminosities calculated from the Ms in Table Asterisks and arrows are defined as in Figure [12j 
We find strong correlations between both C II] and Mg II and L acc . Fits to the data are given in 
Equations [15] and HH 
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Here, we modeled the shock emission as arising in columns carrying a range of energy fluxes 
(or densities) and covering different fractions of the stellar surface, to simulate the highly 
complex geometry of the accretion region. The shock emission is assumed to be a combination 
of high energy flux (J-) columns which peak in the UV and low J- accretion columns which 
peak at redder, optical wavelengths. By including a variety of accretion columns, optical 
excesses and veiling at 1/iin can be explained. Also, the large filling factors of low energy 
columns are in better agreement with analysis based on magnetic field geometries. 

Comparing o u r esti mate s of the accretion prope rties to those calculated in previous analyses 



Valenti et al,l (| 19931 ) and iGullbring et al,l (|1998l ) , we found that our assumptions for Ay and 



intrinsic variability in the optical fluxes accounted for the majority of the differences in our 
M or L acc estimates. Two factors in our technique appear to offset each other; 1) the UV 
flux is decreased by using a WTTS template and 2) fitting 1 fim veiling estimates increases 
the red shock excess. 

Low J- accretion columns which can explain 1 fim veiling in high accretors, may be present 
even when no red excess is observed. The intrinsic stellar emission, mainly the photosphere 
in the red and the active chromosphere in the blue, may hide this accretion component at 
low M. We measured the maximum emission of a cool column that can be hidden by the 
stellar emission and determined the filling factor of the cool column. We found large surface 
coverage, even in sources without veiling at 1 /im. In some cases the upper limit on the red 
excess could be equal to the excess flux in the NUV, doubling M. 

We used simultaneous measurements of emission line luminosities and accretion luminosities 
to calibrate correlations between the two. We found clear trends between the luminosities 
of H/3, Ca II K, Mg II and CII]. No correlation is found between EW(Ha) and M, even 
with quasi-simultaneous observations, however the luminosity of Ha is correlated with the 
accretion properties. C II] may be a sensitive tracer of accretion at low M, as it is observed 
in all the CTTS in our sample, even those with no NUV accretion excess, yet absent from 
the spectra of WTTS. 
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Table 1. WTTS Sample and Properties 



Object 


Ay 


Luminosity 


Radius 


Mass 


SpT 




mag 


(L Q ) 


(Re) 


M Q 




HBC 427 


0.0 


0.8 


1.9 


0.8 


K7 


LkCa 19 


0.0 


1.7 


1.6 


1.3 


K0 


RECX 1 


0.0 


1.0 


1.8 


0.9 


K5 


TWA 7 


0.0 


0.5 


1.8 


0.5 


Ml 



Note. Ay and SpT references: HBC 427 

and LkCa 19 dKenvon fc Hartmannl ll995T): R ECX 1 
(jLuhman fc Steeghsl 120041 ): TWA 7 ()Webb et alJll99sh . 
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Table 2. CTTS Sample and Properties 



Object 


SpT 


Ay 


Luminosity 

(L & ) 


Radius 
(Re) 


Mass 
(M Q ) 


Distance 

(PC) 


ry 


>'Y 


AA Tau 


K7 


1.9 


1.0 


2.1 


0.8 


140 


0.3 


0.2 


BP Tau 


K7 


1.1 


1.0 


2.1 


0.8 


140 


0.7 


0.3 


CS Cha 


K6 


0.3 


1.9 


2.7 


0.9 


160 


0.2 


0.0 


CV Cha 


G9 


1.5 


3.1 


2.0 


1.5 


160 


1.1 


0.6 


DE Tau 


M2 


0.9 


0.8 


2.4 


0.4 


140 


0.6 


0.2 


DK Tau A 


K7 


1.3 


1.6 


2.6 


0.7 


140 


0.5 


0.5 


DM Tau 


Ml 


0.7 


0.2 


1.1 


0.5 


140 


0.7 


0.2 


DN Tau 


MO 


0.9 


1.5 


2.8 


0.6 


140 


0.1 


0.0 


DR Tau 


K5 


1.4 


0.4 


1.1 


0.9 


140 


8.1 


2.0 


FM Tau 


MO 


0.7 


0.1 


0.7 


0.6 


140 


>2 


0.8 


bM Aur 


K7 


0.6 


1.2 


2.6 


0.8 


1 A n 

140 


0.2 


0.0 


HN Tau A 


K5 


1.1 


0.7 


1.5 


1.1 


140 


0.8 


0.5 


IP Tau 


MO 


1.7 


0.7 


1.9 


0.6 


140 


0.2 


0.1 


LkCa 15 


K5 


1.1 


0.8 


1.6 


1.1 


140 


0.2 


0.0 


PDS 66 


Kl 


0.2 


0.9 


1.3 


1.1 


86 


0.2 


0.1 


RECX 11 


K5 


0.0 


0.6 


1.4 


1.0 


97 


0.0 


0.0 


RECX 15 


M3 


0.0 


0.1 


0.9 


0.3 


97 


0.8 


0.4 


RW Aur A 


K3 


0.5 


0.5 


1.1 


0.9 


140 


2.0 


0.9 


TWA 3a 


M3 


0.0 


0.4 


1.8 


0.3 


50 


0.0 


0.0 


TW Hya 


K7 


0.0 


0.3 


1.1 


0.8 


56 


0.3 


0.0 


V836 Tau 


K7 


1.5 


1.0 


2.1 


0.8 


140 


0.0 


0.0 



Note. — Ay and SpT references: AA Tau, BP Tau, DE Tau, DN Tau, DR 
Tau, FM Tau GM Aur, HN Tau A, IP Tau, LkCa 15, RW Aur A a nd V836 Tau 



(IFurlan et al.ll201lh : DK Tau A, DM Tau, CS Cha and C V Cha (jFurlan et al 



2009]) : RECX 11 and R ECX 15 (ILuhman fc Stee ghs 2004); TWA 3a and TW 



Hya (TWebb et al.lll999l ): PDS 66 (iMamaiek et all 12002 ). Dista nce references: 



Taur us Molecular Cl oud (IKenvon et al.1 1 1994 ): Ch amaeleon I (IWhittet et al. 



1997) ; y Chamaeleon (jMamaiek et al.lll999l ): TWA (TWebb et al.lll999l): P DS 66 



(jMamaiek et al.ll2002l ). ry and ry values were taken from Edwards et al.1 (J2006J) 
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and lHartmann fc Kenvonl (|l990l ) except for those sources identified in Section 
3.2i ry and ry from the literature are in regular font, whereas those determined 
from new observations or found assuming a relation between ry and ry are in 
italics. 
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Table 3. Log of Observations 



Object 


RA 


DEC 


Telescope/ Instrument 


Date of Obs 




fJ2000) 


fJ2000) 






AA Tau 


04 34 55.42 


+24 28 52.8 


HST/ STIS G230L/G430L 


2011-01-07 








CTIO/ SMARTS RC Spectrograph 


2010-12-31 








CTIO/ SMARTS RC Spectrograph 


2011-01-02 








CTIO/ SMARTS RC Spectrograph 


2011-01-04 








CTIO/ SMARTS RC Spectrograph 


2011-01-05 


BP Tau 


04 19 15.86 


+29 06 27.2 


HST/ STIS G230L 


2002-01-12 


CS Cha 


11 02 25.20 


-77 33 36.3 


HST/ STIS G230L/G430L 


2011-06-01 








CTIO/ SMARTS RC Spectrograph 


2011-05-30 








CTIO/ SMARTS RC Spectrograph 


2011-05-31 








CTIO/ SMARTS RC Spectrograph 


2011-06-02 








Magellan/ MIKE 


2012-02-16 








VLT/ CRIRES 


2010-12-23 


CV Cha 


11 12 27.65 


-76 44 22.1 


HST/STIS E230M/G430L 


2011-04-13 








CTIO/ SMARTS RC Spectrograph 


2011-04-08 








CTIO/ SMARTS RC Spectrograph 


2011-04-13 








CTIO/ SMARTS RC Spectrograph 


2011-04-15 








CTIO/ SMARTS RC Spectrograph 


2011-04-17 








CTIO/ Bench- Mounted Echelle Spectrograph 


2011-04-12 








Magellan/ MIKE 


2010-03-10 


DE Tau 


04 21 55.69 


+27 55 06.1 


HST/ STIS G230L/G430L 


2010-08-20 








CTIO/ SMARTS RC Spectrograph 


2010-08-15 








CTIO/ SMARTS RC Spectrograph 


2010-08-18 








CTIO/ SMARTS RC Spectrograph 


2010-08-20 








CTIO/ SMARTS RC Spectrograph 


2010-08-21 


DK Tau A 


04 30 44.25 


+26 01 24.5 


HST/ STIS G230L/G430L 


2010-02-04 


DM Tau 


04 33 48.74 


+18 10 09.7 


HST/ STIS G230L/G430L 


2010-08-22 








CTIO/ SMARTS RC Spectrograph 


2010-08-15 








CTIO/ SMARTS RC Spectrograph 


2010-08-18 








CTIO/ SMARTS RC Spectrograph 


2010-08-20 








CTIO/ SMARTS RC Spectrograph 


2010-08-21 








Magellan/ MIKE 


2011-01-04 








VLT/ CRIRES 


2011-08-11 
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Table 3 — Continued 



Object 


RA 


DEC 


Telescope/ Instrument 


Date of Obs 




(J2000) 


(J2000) 






DN Tau 


04 35 27.44 


+24 14 59.1 


HST/ STIS G230L/G430L 


2011-09-10 








CTIO/ SMARTS RC Spectrograph 


2011-09-07 








CTIO/ SMARTS RC Spectrograph 


2011-09-11 


DR Tau 


04 47 06.22 


+16 58 42.6 


HST/ STIS G230L/G430L 


2010-02-15 


FM Tau 


04 14 13.56 


+28 12 48.8 


HST/ STIS G230L/G430L 


2011-09-21 








CTIO/ SMARTS RC Spectrograph 


2011-09-21 








CTIO/ SMARTS RC Spectrograph 


2011-09-27 








VLT/ CRIRES 


2010-12-02 


GM Aur 


04 55 10.98 


+30 21 59.1 


HST/ STIS G230L/G430L 


2010-08-19 


HBC 427 


04 56 02.02 


+30 21 03.2 


HST/ STIS G230L/G430L 


2011-03-30 


HN Tau A 


04 33 39.37 


+17 51 52.1 


HST/ STIS G230L/G430L 


2010-02-10 


IP Tau 


04 24 57.14 


+27 11 56.4 


HST/ STIS G230L/G430L 


2011-03-21 








VLT/ CRIRES 


2010-11-13 


LkCa 15 


04 39 17.73 


+22 21 03.8 


HST/ STIS G230L 


2003-02-13 








Magellan/ MIKE 


2011-01-04 








VLT/ CRIRES 


2010-02-03 








VLT/ CRIRES 


2010-02-05 


LkCa 19 


04 55 36.97 


+30 17 55.0 


HST/ STIS G230L/G430L 


2011-03-31 


PDS 66 


13 22 07.45 


-69 38 12.6 


HST/ STIS G230L/G430L 


2011-05-23 








CTIO/ SMARTS RC Spectrograph 


2011-05-16 








CTIO/ SMARTS RC Spectrograph 


2011-05-20 








CTIO/ SMARTS CHIRON 


2011-05-20 








CTIO/ Bench-Mounted Echelle Spectrograph 


2011-05-22 


RECX 1 


08 36 56.12 


-78 56 45.3 


HST/ STIS G230L/G430L 


2010-01-22 








VLT/ CRIRES 


2011-05-30 


RECX 11 


08 47 01.28 


-78 59 34.1 


HST/ STIS G230L/G430L 


2009-12-12 








CTIO/ SMARTS RC Spectrograph 


2009-11-27 








CTIO/ SMARTS RC Spectrograph 


2009-12-15 








CTIO/ SMARTS RC Spectrograph 


2009-12-19 








Magellan/ MIKE 


2010-03-10 








VLT/ CRIRES 


2011-05-30 


RECX 15 


08 43 18.43 


-79 05 17.7 


HST/ STIS G230L/G430L 


2010-02-05 
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Table 3 — Continued 



Object 


RA 


DEC 


Telescope/ Instrument 


Date of Obs 




(J2000) 


(J2000) 












CTIO/ SMARTS RC Spectrograph 


2010-02-01 








CTIO/ SMARTS RC Spectrograph 


2010-02-04 








CTIO/ SMARTS RC Spectrograph 


2010-02-06 








Magellan/ MIKE 


2010-03-11 


RW Aur A 


05 07 49.51 


+30 24 04.8 


HST/ STIS G230L/G430L 


2011-03-25 


TWA 3a 


11 10 27.80 


-37 31 51.2 


HST/ STIS G230L/G430L 


2011-03-26 








CTIO/ SMARTS RC Spectrograph 


2011-03-18 








CTIO/ SMARTS RC Spectrograph 


2011-03-22 








CTIO/ SMARTS RC Spectrograph 


2011-03-25 








CTIO/ SMARTS RC Spectrograph 


2011-03-29 








Magellan/ MIKE 


2010-03-10 


TWA 7 


10 42 29.94 


-33 40 16.7 


HST/ STIS G230L/G430L 


2011-05-05 


TW Hya 


11 01 51.95 


-34 42 17.7 


HST/ STIS G230 


2002-05-10 


V836 Tau 


05 03 06.62 


+25 23 19.6 


HST/ STIS G230L/G430L 


2011-02-05 








CTIO/ SMARTS RC Spectrograph 


2011-01-28 








CTIO/ SMARTS RC Spectrograph 


2011-02-01 








CTIO/ SMARTS RC Spectrograph 


2011-02-03 








CTIO/ SMARTS RC Spectrograph 


2011-02-08 
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Table 4. Literature Ays for Taurus Sources 



Object 


F09/11 


KH95 


G98 


V93 


Fll 


M13 


AA Tau 


1.9 


0.5 


0.7 


1.3 


1.3 




BP Tau 


1.1 


0.5 


0.5 


0.9 


1.8 


0.6 


DE Tau 


0.9 


0.6 


0.6 


1.7 




0.9 


DK Tau A 


1.3 


0.8 


1.4 


1.2 


1.8 




DM Tau 


0.7 


0.0 




0.1 






DN Tau 


0.9 


0.5 


0.3 


0.5 






DR Tau 


1.4 






1.0 


1.8 


2.0 


FM Tau 


0.7 


0.7 




0.8 






GM Aur 


0.6 


0.1 


0.3 


0.5 






HN Tau A 


1.1 


0.5 


0.7 


1.0 


3.1 




IP Tau 


1.7 


0.2 


0.3 








LkCa 15 


1.1 


0.6 










RW Aur A 


0.5 






1.2 






V836 Tau 


1.5 


0.6 








1.4 



Note. - F09/11 


(Furlan et a 


2009. 


2011 


), KH95 


(Kenvon & Hartmann 


1995), G98J( 


Uullbrine et al. 


1998V 



V93 (Valcntietal 



(McClure et al 



1 22^ FIJ (Fischor H al 



in prepl ) 



2011 



M13 
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Table 5. Results from Multi-Component Model Fits to UV and Optical Spectra 



Object 


/(io 10 ) 


/(3 x 10 10 ) 


/(io 11 ) 


/(3 x 10 11 ) 


/(io 12 ) 


f total 


ry 


M (Mo yr^ 1 ) 


A A m 

AA Tau 














0.002 


0.002 


0.0 


1.5x10 


BP Tau 








0.02 





0.002 


0.022 


0.2 


2.9xl0 -8 


CS Cha 


0.008 





0.003 








0.011 


0.0 


5.3xl0~ 9 


CV Cna 





0.4 





0.02 





0.42 


0.6 


5.9xl0 -8 


DE Tau 


0.1 





0.0008 


0.0007 





0.1 


0.3 


2.8x10 


DK Tau A 


0.01 





0.005 


0.005 





0.02 


0.1 


3.4x10 


DM Tau 


0.08 








0.001 


0.0006 


0.082 


0.3 


2.9xl0~ 9 


DN Tau 











0.002 





0.002 


0.0 


l.OxlO" 8 


DR Tau 








0.3 


0.07 


0.006 


0.37 


1.6 


5.2xl0~ 8 


FM Tau 





0.07 








0.001 


0.071 


0.8 


1.2xl0~ 9 


GM Aur 








0.001 


0.003 





0.004 


0.0 


9.6xl0~ 9 


HN Tau A 











0.01 


0.004 


0.014 


0.1 


1.4xl0~ 8 


IP Tau 














0.001 


0.001 


0.0 


7.2xl0~ 9 


LkCa 15 








0.01 


0.0007 


0.0001 


0.011 


0.1 


3.1xl0~ 9 


PDS 66 














0.0001 


0.0001 


0.0 


1.3xl0~ 10 


RECX 11 








0.001 








0.001 


0.0 


1.7X10" 10 


RECX 15 


0.02 


0.001 





0.0007 


8xl0~ 5 


0.022 


0.3 


8.0xl0" 10 


RW Aur A 





0.03 


0.2 








0.23 


0.8 


2.0xl0~ 8 


TWA 3a 














8xl0~ 6 


8xl0~ 6 


0.0 


9.8X10" 11 


TW Hya 











0.0013 


0.0013 


0.0026 


0.0 


1.8xl0~ 9 


V836 Tau 














0.0002 


0.0002 


0.0 


l.lxl0~ 9 



Note. — The values in parentheses represents the energy flux (J 7 ) of each column in erg s 1 cm 
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Table 6. Accretion Rates from Literature for Taurus Objects 



Object 


M(G98) a 
(M yr" 1 ) 


M(V93) 
(M yr" 


b 

') 


AA Tau 


3.3x10" 


-9 


7.1x10" 


-9 


BP Tau 


2.9x10" 


-8 


2.4x10" 


-8 


DE Tau 


2.6x10" 


-8 


1.8x10" 


-7 


DK Tau A 


3.8x10" 


-8 


6.1x10" 


-9 


DM Tau 






2.9x10" 


-9 


DN Tau 


3.5x10" 


-9 


1.3x10" 


-9 


FM Tau 






7.7x10" 


-9 


GM Aur 


9.6x10" 


-9 


7.4x10" 


-9 


HN Tau A 


1.3x10" 


-9 


3.9x10" 


-9 


IP Tau 


8.0x10- 


10 






RW Aur A 






3.3x10" 


-7 



Note. — 

a Gullbring et"all d 19981 1 
MValenti et al.1 (|l993h 
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Table 7. Maximum / and M Assuming Hidden Accretion Emission 



Object 


M a (M© yr- 1 ) 


r 


M b (M yr- 1 ) 


f 


AA Tau 


1.5xl(T 8 


0.002 


1.6xl0~ 8 


0.08 


BP Tau 


2.9xl(T 8 


0.022 


3.3xl0~ 8 


0.2 


DN Tau 


l.OxKT 8 


0.002 


1.7xl0~ 8 


0.06 


DR Tau 


5.2xl(T 8 


0.37 


5.6xl0~ 8 


0.8 


GM Aur 


9.6xl(T 9 


0.004 


1.3xl0~ 8 


0.04 


HN Tau A 


1.4xl0" 8 


0.014 


1.5xl0" 8 


0.06 


IP Tau 


7.2xl0~ 9 


0.001 


9.4xl0~ 9 


0.03 


LkCa 15 


3.1xl0~ 9 


0.011 


3.6xl0~ 9 


0.03 


PDS 66 


1.3xl(T 10 


0.0001 


3.8xl0~ 10 


0.02 


RECX 11 


1.7xl0- 10 


0.001 


2.6xl0- 10 


0.006 


TWA 3a 


9.8xl0~ n 


8 x 10~ 6 


3.4xl0~ 10 


0.002 


TW Hya 


1.8xl0~ 9 


0.0026 


2.3xl0~ 9 


0.05 


V836 Tau 


l.lxl0~ 9 


0.0002 


2.6xl0~ 9 


0.02 



Note. - 

a Values which give the best Xred ^ of the model to the UV and 
optical data (see Table [5]). 

b Values calculated when allowing for a low T column which pro- 
duces emission not detectable above the intrinsic stellar emission. 



